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A cluster of three spins with single-axis anisotropic exchange coupling exhibits a range of classical 
behaviors, ranging from regular motion at low and high energies to chaotic motion at intermediate 
energies. A change of variable makes it possible to isolate total angular momentum around the z-axis 
(a conserved quantity) and it's associated cyclic variable from two non-trivial degrees of freedom. 
This clarifies the interpretation of Poincare sections, and causes permutation symmetries of the 
system to manifest as rotation symmetries in the new coordinates. 

The three-spin system has four families of periodic orbits (multiple instances of each orbit exist 
because of permutation symmetry.) Analysis of spin waves predicts their periods in the low energy 
(antiferromagnetic) and high energy (ferromagnetic) limits and also can be used to determine the 
stability properties of certain orbits at intermediate E = ~1. The system also undergoes interesting 
changes in the topology of the energy surface along particular curves in energy-parameter space, for 
instance, when two pieces of the energy surface surrounding the two antiferromagnetic fixed points 
coalesce. Poincaresections produced with a 3-d graphic technique (described in the Appendix) 
illustrate the symmetries of the system and illustrate the transition to chaos at low energies. 

The three spin system turns out to have similarities with the Anisotropic Kepler Problem (AKP) 
and the Henon-Heiles Hamiltonian. An appendix discusses numerical integration techniques for spin 
systems. The quantum manifestations of the structures found in this paper are discussed in [P. A. 
Houle, N. G. Zhang, C. L. Henley, Phys. Rev. B 60 15179 (1999)]. 

PACS numbers: 03.65.Sq, 05.45.-|-b 



I. INTRODUCTION 

This paper is a study of the classical mechanics of 
a cluster of three spins with single-axis anisotropic ex- 
change coupling; the Hamiltonian of the system is 



H = J 



l<i<j<3 



(1.1) 



We can set J — 1 without losing generality (even the 
s ign is arbitrary for our study of the dynamics). @@ Eq. 
(1.1) is a model for spin clusters in triangular antifer- 
romagnets such as NaTi02 and RbFeCla. It was intro- 
duced in Q and is further reported on in ^, ||]. In 
these works, Nakamura studied the level statistics of the 
system, and compared the behavior of energy level fluc- 
tuations as a function of a in regions where the classical 
dynamics was predominantly regular and chaotic. In a 
previous work, we've discussed how the classical struc- 
tures of the system are manifested in its quantum spec- 
trum 0. 

This system is of intrinsic interest as a comparatively 
simple dynamical system possessing high symmetry. We 
were specifically motivated by semiclassics, i.e. the rela- 
tion between classical dynamics and the eigenstates of a 
quantum system. 

The quantum manifestations of classical cha os were in- 
vestigated not only on the model of Eq. (|l.l|), but also 
in a cluster of two interacting spins j9[ |l^, |ll|, |12|, p^ , Q . 
The method of Bohr-Sommerfeld quantization has also 
been applied, but to single spins |^ . 

In recent years, a variety of magnetic molecules 



have been synthesized containing clusters of interacting 
spins 1^. Although commonly approximated as a single 
moment, they do have internal excitations which can- 
not be found exactly by exact diagonalization (since the 
Hilbert space is too large). If the spin cluster divides 
into subclusters, each of which has a moderately long 
moment, these excitations may best be grasped semi- 
classically. The methods of the present paper would be 
a natural starting point for such studies, at least for sys- 
tems with high symmetry. 

Eq. ( pTl] ) conserves total Sz = Yl^i=i^iz- A key as- 
pect of our approach is a change of variable which sep- 
arates the system into two subsystems: (i) A noninte- 
grable two-degree of freedom subsystem, and (ii) a sin- 
gle degree of freedom subsystem for which total Sz is a 
conserved momentum variable. Our change of variable 
makes the Poincare sections comprehensible [in contrast 
to prior work that did not separate subsystems (i) 
and (ii)], and makes it easy to find the fundamental pe- 
riodic orbits of the system. As in previous work on Eq. 
(1.1), we restrict our consideration to the special case 
that Sz — 0; as this paper considers only the classical 
mechanics of the system, we set |S| = 1 without losing 
generality. 

An outline of this paper follows: Section |l| introduces 
our change of variable and shows how the topology of 
the phase space and the symmetries of t he system ap- 
pear in the new coordinates. Section III is a discussion 



of the fixed points, invariant manifolds, and spin waves 
of the three-spin system. Section IV is about changes 



in the topology of the energy surface that occur along 
certain curves in the E — a parameter space. Section ^ 
enumerates the fundamental periodic orbits of the sys- 
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tem and maps the global d yna mics of the system with 
Poincare sections. Section ^ concludes by discussing 
similarities between the the three-spin cluster and the 
well-studied Anisotropic Kepler Problem and the Henon- 
Heiles Hamiltonian. 



II. HEXAGONAL PHASE SPACE 

With three degrees of freedom, the dynamics of an ar- 
bitrary dynamical system is difficult to visualize and un- 
derstand. Fortunately, the three-spin cluster conserves 
total angular momentum around the z-axis, so we can 
separate the problem into two halves: (i) an autonomous 
nonintegrable two degree of freedom Hamiltonian, and 
(ii) an trivial single degree of freedom system which is 
driven by system (i) . Part (i) can be studied in isolation 
from part (ii), however, since (ii) is driven by (i) we must 
solve the motion of (i) before we can solve the motion of 
(ii). 



Section P A| presents a change of variable that sepa- 
rates part (i) from part (ii) - which is necessary to draw 
useful Poincare sections. Next, Section [I B shows how 
the discrete symmetries of the s yste m appear in the new 
coordinates. Then, in Section II C we'll show that the 
27r periodicity of the spin longitudes (0^) manifests as 
hexagonal tiling in the {^a,^b) coordinates. 



A. A change of coordinates 

The conventional canonical coordinates for spin are 
{(j)i,Zi) where is longitude and Zi — Siz- ||l5|] Spin 
space includes both position and momentum and is the 
complete phase space of a spin system - the position 
of the spin vectors at a moment in time completely 
describes the system. The mapping between {(pi, Zi) 
and spin vector components is Si — {Six, Siy, Siz} — 



{(l-z2)icos0„(l-zf) 



2 sm( 



J- 



To isolate total angular momentum around the z axis, 
we make the orthogonal linear transformation 



$0 

where 




(2.1) 



M = — 




\/2 V2 V2 

2 -1 -1 

V3 -Vs 



(2.2) 



Because $o is a cyclic coordinate (does not appear 
in Eq. (2.3)), Zq is conserved: Zq is proportional to 
total spin around the z-axis and $o measures the col- 
lective precession of the three spins around the z-axis. 
We call the two degree-of-freedom system consisting of 
{'^ A,^ B ,Z A,Z b) the reduced system and the complete 



system with three degrees of freedom the full system. 
As the equations of motion for the reduced system do 
not depend on $0, the evolution of $0 can be ignored 
when we study the reduced system; enters only as a 
constant parameter of the reduced system. 

Zq = when total 5^ = is zero (as is always the case 
in this paper), and the Hamiltonian becomes 



H 



(l-a) 



2 +Z|)+/+/_cos(\/2$b) 




in the new coordinates with 



/o 



2 „ 
1 - 3^1, 



f± = \Ii-^[Za±V3Zi 



(2.3) 

(2.4) 
(2.5) 



If we're interested in the evolution of $0, we can study 
the full system by first finding a trajectory of the reduced 
system and then solving the remaining time-dependent 
equation of motion for $0. A periodic orbit of the re- 
duced system may or may not be a periodic orbit of the 
full system; a periodic orbit of the reduced system is a 
periodic orbit of the full system only if $0 changes by 
an integer multiple of 27r/\/3 per orbit of the reduced 
system. 



B. Symmetries in $, Z coordinates 

Structures in phase space, such as invariant manifolds 
and fundamental periodic orbits, reflect symmetries of 
a dynamical system. Therefore, it's essential to under- 
stand the discrete symmetries of a system in order to 
characterise structures in phase space. 

Permutations of the identities of the spins are an im- 
portant set of discrete symmetry operations for the three 
spin cluster. For instance, if we swap S2 and S3 the 
Hamiltonian Eq. (1.1) is unchanged. The three spin 
cluster is the = 3 case of a simplex, in which all of the 



spins and the bonds between the spins are interchange- 
able. The associated symmetry group is the permutation 
group P3 which requires two generators. The first gen- 
erator is T, the translation operator. T maps Si — s- S2, 
S2 — > S3, and S3 — > Si. Viewed in ($, Z) coordinates, T 
is a rotation of both the {^a,^b) and {Za,Zb) planes 
by 27r/3 or 



(2.6) 
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with 




(2.7) 



the same transformation over Z. The exchange operators 
are a subset of i-3: E^, where i is a spin index, exchanges 
the other two spins. Ei is 




with 



El 




(2.8) 



(2.9) 



in ($, Z) coordinates. (The same transformation matrix 
also acts on the Z coordinates). Any Ei can be chosen 
for a second generator - T and Ei are a complete set of 
generators for P3. 



C. The tiling of $, Z coordinates 

Unlike the phase space of particle systems, which is 
infinite in extent, spin space is compact and periodic in 
0i. Therefore, spin trajectories exist which have no ana- 
log in a particle system. For instance, a spin can precess 
around the z-axis and come back to its initial position 
without the sign of ever changing. Also, a spin can 
pass directly over the north pole [zi — > 1) at which point 

jumps discontinuously by tx. 

The transformation Eq. (2J^) also changes the ap- 
pearance of the boundaries and connectivity of the 
phase space; the periodicity of the coordinates causes 
{^A^ ^b) to be periodic on a hexagonal lattice while the 
condition \zi\ < 1 restricts (Za,Zb) to the interior of a 
hexagon. 

Because the the coordinates are 2tt periodic the 
transformation (pi ^ (pi + 2tt does not change the state 
of the system. This periodicity looks different in the 
(<i>^, $b) plane: the transformation —^ (f)i + 27r maps 
to to a translation in {^a, ^b) space of (^j 0). Adding 

2tt to (p2 maps to a translation of ■2^(— 1, -^/s) and adding 
27r to 03 maps to a translation of ■2i(— 1, — \/3) - related 
by the operator T, these three vectors form the lattice 
vectors of a hexagonal lattice, see Fig. |l|. 

Since the z's and Z's are related by the same linear 
transformation that relates the 0's and <i>'s, the domain 
of valid {Za,Zb) is hexagonal. One boundary of the 
hexagon is where spin 1 is at the north pole, zi = I, 
since Zi € [—1,1]. At that point, Za = When spin 1 

is at the south pole, zi = — 1 and Za = The rest 

of the boundaries can be found by rotating the zi = ±1 
boundaries by ±120° in the {Za, Zb) plane. 



A remaining detail is how the Z trajectory appears 
as a spin passes through a pole. For example, when spin 1 
hits the north pole, the the trajectory strikes the bound- 
ary of the {Za,Zb) hexagon a,t Za = Although 
there is no discontinuity in the z or Z coordinates, the 
01 coordinate jumps discontinuously by tt. As seen in 
{^a,^b) coordinates, the coordinate ^a jumps by 
As z'l changes discontinuously at this point, the projec- 
tion of the trajectory seems to "bounce" off the boundary 
in {Za, Zb-) 



III. FIXED POINTS, INVARIANT MANIFOLDS 
AND SPIN WAVES 

A. Fixed Points 

Certain phenomena of the three spin system, such as 
fixed points, invariant manifolds and spin waves, can 
be studied without numerical integration. Section [II A 



concerns the fixed points of the three-spin system in 
whic h all spins lie in the equatorial plane. Next, Sec- 
tion 



IIIB describes invariant manifolds of the system 



two-dimensional subspaces on which the dynamics are 
reduced to a single degree of freedom. Finally, Section 
III C develops a linear expansion around the fixed points 



found in Section III A to derive the frequency of spin 



wave excitations in their phase-space vicinity. 

Fixed points are points in phase space where $j = 
Zj = for j G {0, A, B} There are three families of fixed 
points of the system in which the spins lie in the equa- 
torial plane {Za = Zb = 0.) The locations of these fixed 
points in {^a, ^b) space are plotted on Fig. 0. The fer- 
romagnetic state {FM) is the global energy maximum 
{E = 3) with all three spins pointing together. The FM 
state lies at (0,0) in the {^a,^b) plane. The two an- 
tiferromagnetic states {AFMj^, and AFMl) are located 
at (0, ±^^^) in {^a,^b) plane and are mirror images 
of each other. The antiferromagnetic states are global 
energy minima {E = —1.5) with the spins splayed 120° 
apart. There are also three antiparallel configurations 
{Ai where z is a spin index i S {1,2,3}) where two spins 
are coaligned while the other spin (spin i) points in the 
opposite direction; here E = —1. Ai lies at (0, V^n) while 
A2 and A3 lie at {± y/3/2Tr , ^/T/2t:) in the {^a,<^b) 
plane - the operator T transforms into another by the 
operator T, a 120° rotation in the {^a,^b) plane. The 
F, AFM and A fixed points are fixed points of both the 
reduced and full systems. 



B. Invariant Manifolds 

If a subspace of the phase space is an invariant mani- 
fold, the time evolution of the system will remain in that 
subspace if its initial state lies in that subspace. As the 
invariant manifolds of the reduced system are two dimen- 
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sional, dynamics on the invariant manifolds possess only 
a single degree of freedom. Therefore, a family of peri- 
odic orbits lives on each invariant manifold, which are 
discussed further in Section The three spin system 
has two kinds of invariant manifold: the stationary spin 
manifolds and the counterbalanced manifolds. For either 
kind of manifold, the motion of one spin is different from 
the other two; with the operator T one can find three 
manifolds of each type, related to one another by a 120° 
rotation in the {^a,^b) or {Za,Zb) planes. 

The = 'Z'A = subspace is one stationary spin 
manifold. On this manifold, spin 1 lies in the equatorial 
plane and remains stationary while the other two spins 
execute roughly circular motions in opposite directions. 
Motion on the stationary spin manifolds can be modelled 
with a one-spin system: pointing spin 1 along the x-axis, 
the constraint = Za = combined with ■^'z = 
implies that Ssx = S2x, S^y = -S2y and 6*32 = -S^iz- In 
spin vector form, the reduced Hamiltonian is 



Hss = 25, + Sl- Si - (1 - a)Sl 



(3.1) 



where 8 = 82. 

The other family of invariant manifolds are the coun- 
terbalanced manifolds in which two spins move together 
in a direction opposite to the other spin; the counter- 
balanced manifold with spin one the special spin is the 
subspace $b = Q,Zb — 0. With the arbitrary choice 
of $0 = 0, the following constraints apply: 82 = 83, 
S\y = — 25*2^ and Siz — —25*22. The remaining con- 
straint, on 523;, is determined by the total length con- 
straint I 81 1= 1, which imphes S2X = (+^3 + <5'L) /2- 
With 8 = Si the Hamiltonian reduces to 



H = J 



1 + Sx^^+Sl 



Si - St 
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(3.2) 



In the case of tr = 0.5, for which the quantum mechanics 
have been extensively studied (see [Q), the analytic value 
of Tfm = 2.294 agrees with the limit of the periods of all 
fundamental orbits (computed by numerical integration, 
see 0) as 3. 

A similar expansion is possible around either AFM 

ground state , where ^a and 8b = $_b — are small. 
We obtain 



Ha 



FM 



Ea 



FM 



51) 



1 



zl) 



(3.5) 



Here there are two degenerate spin waves with period 



Tafm{(^) 



2 27r 
3 



(3.6) 



Approaching the isotropic case, ti — s- 0, Tafm becomes 
infinite. This coincides with the exact solution for a — Q 
in which all three spins precess around the total spin 
vector Q at a rate proportional to the length of the to- 
tal spin vector: at AFM , the total spin and the rate of 
spin precession are both zero. When a = 0.5, the peri- 
ods of all fundamental orbits converge to Tafm = 7.255 
as -E — > — 1.5 as predicted by Eq. ( ^.6[ ). As the spin 
wave frequency drops to zero, the zero point energy of 
the quantum ground state also drops to zero, converging 
on the classical ground state energy as is observed in [0 . 

Although Ai is a saddle point rather than a ground 
state, it is still possible to linearize the Hamiltonian in its 
vicinity. Let ^b and 5 A be small, where 5 A — ^a — ^A, 
and 3>A = -v/2737r. Then, 



Ha ~ Ea 



6 + 2^^+2^- 



(1-a) 



Zi 



$1 



(3.7) 



C. Spin Waves 

Small-energy excitations of a spin system understood 
in terms of the linearized dynamics around a fixed point 
are spin waves. This section is a study of the linearized 
dynamics around the ferromagnetic {FM), antiferromag- 
netic {AFM) and antiparallel fixed points {A). The term 
'spin wave' is usually used to refer to excitations of a 
ground state, but the concept remains useful at points 
such as the antiparallel fixed points which are saddles of 
the energy function. 

The linearization of Eq. ( ^.3D near the the ferromag- 
netic (FM) fixed point is 



H 



FM 



3-cr 



[Zl 



zl) 



There are two degenerate spin waves with period 

27r 



l-FM 



{a) = 



^3(3 - <J) 



(3.3) 



(3.4) 



The first set of terms in Eq. (3.7) depends on Za and 
^A and the second set depends on Zb and $3. The first 
set in Eq. (3.7) describes positive energy spin waves that 
live on the Zb = =0 counterbalanced manifold with 
period 



27r 



V1 + 3(T 



(3.8) 



which is Ta.c = 3.975 when a — 0.5. The second set 
describes negative energy spin waves that live on the 
Za = ^A = stationary spin manifold with period 



Ta,ss{<7) 



2n 



(3.9) 



which is Ta,ss — 8.885 when a = 0.5. Corners of the 
counterbalanced and stationary spin manifolds touch at 
right angles at A. We will later take advantage of this 
to compute the stability properties of the stationary spin 
and counterbalanced orbits in Section 
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IV. THE TOPOLOGY OF THE ENERGY 
SURFACE 

Unlike the canonical cases of two-degree of freedom 
Hamiltonian dynamics, such as the Anisotropic Ke- 
pler Problem (AKP) and the Henon-Heiles Hamiltonian, 
the three-spin system exhibits nontrivial changes in the 
topology of the energy surface at certain energies. There 
are three transition energies: (i) the coalescence energy 
Ec{(j)\ (ii) the antiparallel energy E a = — 1; and (iii) the 
polar energy Ep{<T). Fig. ^depicts the transition energies 
as a function of a while Fig. ^ illustrates the topology 
changes for cr = 0.5. As we increase energy from the 
ground state E — —1.5, transition (i) always occurs first. 
If < cr < 2/3, transition (iii) occurs before transition 
(ii), otherwise when 2/3 < a < 1 (ii) occurs before (iii). 



A. The coalescence transition 

At low energies (near the AFM fixed points) an en- 
ergy barrier separates the two antiferromagnetic ground 
states. Therefore, the energy surface is composed 
of two disconnected parts. Those parts become con- 
nected when E — Ec{a). Using polar coordinates, S = 
(sin0cos(/), sin sine/), cos 0), the surfaces first touch at 
the saddle points on the = line on Fig. || (one is 
hidden behind the sphere.) Thus, Ecicr) is found by con- 
sidering Eq. (3.1), the single-spin Hamiltonian for the 
stationary-spin manifold. The saddle lies on the <f> — Q 
line, at the point where dH/d9 = 0, or 



sin 6 



1 



2-cr 



(4.1) 



Substituting this back into Eq. (|3.lD, the saddle-point 
energy is 



3 - 3(7 -h cr^ 
cr-2 



(4.2) 



The coalescence occurs between (a) and (b) in Fig. 
^. An interesting quantum manfestation of the coales- 
cence transition was observed in ||l[, the tunnel splitting 
of quantum levels as Ec is approached from below. 



B. The antiparallel transition 

At Ec, a set of necks come into existence that connect 
the two lobes of the energy surface that were discon- 
nected at energies below Ec (See Fig. IJb) when E = —1 
these necks fuse, changing the connectedness of the en- 
ergy surface again (See Fig. ^.) This is the antiparallel 
trans ition - the antiparallel fixed point Ai (see Section 
III A ) lies in the center of the large hole in Fig. Two 
families of periodic orbits disappear at this transition, 
including one branch of stationary spin orbits approach- 
ing from ii' < — 1 as well as the counterbalanced orbits 



approaching from E > ~1. (One aspect of Fig. |jb is 
deceptive. Being a three-dimensional cut out of a four- 
dimensioanl space, it fails to show two other pairs of con- 
necting necks that surround the A2 and A3 antiparallel 
fixed points - for a total of six necks.) 



C. The polar transition 

The system attains extreme energies (both minimum 
and maximum) only when the spins lie in the equato- 
rial plane. As a result, there is both a minimum and a 
maximum energy at which one spin can point at a pole, 
which is a saddle point in the full phase space. These 
are the upper and lower polar transition energies - these 
thresholds are found by pointing one spin, say spin 1, 
at the north pole and finding the maximum and mini- 
mum energy configurations. Setting Za = -\/3/2, in Eq. 
( ^3pl| ) we get 



H = - 



2 
1 

+ 4 



3 , 72 



20Z| 



AZ% cos (2$b) . 



(4.3) 



which has a maximum at H„ 



7(7 and a minimum 



at Hmin = fcr ~ The polar transition occurs between 



panels (c) and (d) in Fig. IJ: at this point the energy 
surfaces touch the enclosing hexagonal prism, forming a 
network of necks connecting the energy surface to itself. 



D. The classical density of states 

Changes in the topology of the surface section have an 
interesting effect on the classical and quantum densities 
of states. The weighted area of the energy surface, 

PciE) = J d^Ad^BdZAdZeSiE - H{H, E)) (4.4) 

is the classical density of states, since it is proportional 
to the quantum density of states. |l^ Fig. ^ is a plot 
of the classical density of states for the reduced system 
as a function in energy. We observe two interesting fea- 
tures: first, a discontinuity in the slope of ac{E) at the 
coalescence transition (This is (a) in Fig. ^.) Second, 
the density of states is apparently flat between the lower 
polar transition Ep{(j) and the antiparallel transition - 
although we don't have an analytic understanding of the 
flat spot, numerical evidence suggests that it is exactly 
flat. 



V. FUNDAMENTAL PERIODIC ORBITS AND 
GLOBAL DYNAMICS 

This section presents the main results we've deter- 
mined from numerical integration of the equations of 
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motion: a map of the fundamental periodic orbits of 
the three-spin system for and Poincare sections depict- 
ing the global dynamics of the system for a — 0.5. Like 
any chaotic system, the three spin system has an infinite 
number of periodic orbits. However, a few short period 
orbits form the skeleton of the system's dynamics. Four 
of these are known; the stationary- spin orbit, the coun- 
terbalanced orbit, the three-phase and the unbalanced or- 
bit. Fig. 1^ plots the energy-time curves of the four orbits 
for a = 0.5. (Spin trajectories for the four orbit types are 
visualized in Fig. 1 of 



Sections 



the stationary-spin manifold, the period of the stationary 
spin orbits goes to infinity as E Ec from either side, 
with an observable effect on the quantum mechanical or- 
bit spectrum. 

In the case of cr = 0.5, the stationary spin orbit is 
unstable for E > —1. For —l<E<Ec, the inner branch 
of the stationary spin orbit is stable and the outer branch 
is unstable. The left and right branches are stable as — > 
— 1.5 but become unstable as the energy increases and 
chaos becomes widespread (the orbit does momentarily 
regain its stability near E = —1.22.) Fig. |^ is a plot 
V A| - |VC| discuss the stationary-spin, coun- of the stability parameter p = A + A^^ versus energy 

for E < —1, where A and A^^ are eigenvalues of the 
transverse stability matrix; |p| < 2 for a stable orbit and 
and IpI > 2 an unstable orbit. Il^] 



terbalanced, thr ee-ph ase and unbalanced orbits respec- 
tively. Section VE and |VF| discuss global dynamics 
near the ferromagnetic {E 3 ) and antifcrromagnetic 
(E —1.5) ends. Section V G points out how symme- 



tries of the Hamiltonian manifest in its classical dynam- 



B. Counterbalanced orbits 



A. Stationary spin orbits 

The stationary spin orbits are simple to study because 
any poi nt on a stationary spin invariant manifold (see 
Section III B ) lies on a stationary spin orbit. As there are 
three stationary spin manifolds, there are three families 
of stationary spin orbits related by symmetry. 

For a stationary spin orbit, one spin (say, spin 1) is 
stationary in the equatorial plane, while the other two 
spins move in distorted circles, 180° out of phase. Fig. 
|5| is a plot of the trajectories of one of the moving spins, 
based on the single-spin Hamiltonian Eq. (3.1). In the 
range —1<E<3, each family of stationary spin orbits 
has a single branch; trajectories on the single-spin sphere 
are concentric distorted circles centered around the FM 
fixed point. Between E = —I and E = Ec, two branches 
of periodic orbits exist: the outer branch, still centered 
around the FM fixed point, and the inner branch, cen- 
tered around the A fixed point. Below the coalescence 
energy Ec{a), the orbits reorganize into a different pair 
of branches (left and right), one centered around each 
antifcrromagnetic ground state. 

The stationary spin orbit runs along the Za = seam 
on the slice of the energy surface visualized in Fig. ^. 
Fig. ^d and c represent the case where E > —1 and only 
one branch of the orbit exists. In Fig. the outer 
branch runs along the outside of the surface while the 
inner branch rungs along the inside of the hole in the 
surface. Finally in Fig. left and right branches of the 
the stationary spin orbit exist on two separate lobes of 
the energy surface. 

The energy-time curve for the stationary spin orbits 
is seen in Fig. 0. Although one would expect that the 
periods of the left and right branches in the E < Ec 
regime are the same (because they are related by reflec- 
tion symmetry,) it's a bit surprising that the period T{E) 
of the inner and outer branches in the E^ < E < —1 is 
also the same. Since the coalescence separatrix intersects 



The counterbalanced orbit family is also easy to study 
because, like the stationary spin family, it lives on an in- 
variant manifold. Counterbalanced orbits exist only in 
the range E > —1, and at least for < a < 2, the 
counterbalanced orbit is always stable. As is the case for 
counterbalanced manifolds, there are three counterbal- 
anced orbits related by the symmetry T. 



The spin wave analysis of Section III C can be applied 



to the stability properties of the inner branch stationary 
spin and counterbalanced orbits in the limit E ~^ I. A 
counterbalanced orbit (approaching from E > —1) and 
the inner branch of a stationary spin orbit touch at each 
antiparallel fixed point. Using the linearization around 
the antiparallel fixed point Eq. ( |3.7D , we can establish 
that both orbits are stable, and compute the limiting 
value of the stability exponent p dX E ^ —1 for both 
orbits. 

Because there are two distinct frequencies in the lin- 
earized dynamics around the antiparallel fixed point, Eq. 
(|3.7|) , dynamics in the vicinity of the antiparallel fixed 
point are structurally stable and, for close enough ener- 
gies, should be similar to the linearized behavior. Focus- 
ing attention on the counterbalanced orbit, the degree of 
freedom orthogonal to the counterbalanced orbit is the 
stationary spin orbit - therefore the counterbalanced or- 
bit is stable as E —1. In one circuit of the counter- 
balanced spin orbit, a slightly displaced trajectory winds 
around the orbit at the frequency of the stationary spin 
spin wave. The winding number is the ratio of the peri- 
ods of the two orbits, or 



Tea 



1 



1 + 3cr 



For a = 



(5.1) 



2 cos 61 



2 , = . The stability parameter p 
equals —1.891, and agrees with the result obtained by 
numerical integration (see Fig. |^). Repeating this anal- 
ysis for the stationary spin orbit, we obtain p — > 0.175 in 
agreement with Fig. 0. 
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C. Three phase orbits 

The three phase orbits do not he on an invariant mani- 
fold and thus have richer behavior than the previous two 
famihes of orbits. Unhke the stationary spin and counter- 
balanced orbits for which <$o>= 0, three phase orbits 
can exhibit precession, a secular trend in $0 a-nd there- 
fore can be periodic orbits of the reduced system but not 
the full system. The three phase orbit undergoes a pitch- 
fork birfucation at transition energy Eb (Eh « —0.75 for 
a = 0.5.) Above Ei,, a single branch of non-precessing 
orbits exists, but below the bifurcation three branches of 
three phase orbits exist: a non-precessing unstable orbit 
and two stable orbits for which <<i>o>^ with opposite 
signs. This transition is visible on curve (c) of Fig. 0. 
The unstable orbit exists for a small energy below Ei, but 
soon disappears when the spin trajectory intercepts the 
poles of the spin spheres. 

The three-phase orbits are so called because, in three- 
phase orbits, the spins each execute an identical circuit 
around a distorted circle, each 120° out of phase - much 
like the currents used in three-phase AC power transmis- 
sion. The multiplicity of the three-phase orbit is different 
from the previous two: Above Eb there are two three- 
phase orbits, one in which the spins rotate clockwise and 
another with counterclockwise rotation. Below Eb and 
the demise of the nonprecessing orbit, there are a total 
of four: each AFM ground state has it's own pair, one 
member of which has <$o> positive and the other <$o> 
negative. Fig. |l^ is a plot of the per-orbit precession rate 
of the three-phase orbit below the pitchfork bifurcation; 
note that the precession rate converges on 2tt (effectively 
zero) a.s E ^ —1.5. 

The three phase family is more difficult to study than 
the previous two, because we must search the energy sur- 
face for i t. It 's stil l qu ite straightforward, for as seen in 
Sections VE and VF the three spin orbit lies on the 
Zb = line of the — surface of section and can be 
found by a one-dimensional search. The stability expo- 
nent p of the three-spin orbit can be seen in Fig. H. 



D. Unbalanced Orbits 

The unbalanced orbits, a family of unstable periodic 
orbits, exist between E = —1.5 and E — Ep. Unbal- 
anced orbits do not lie on a symmetric manifold and do 
not precess. The unbalanced family corresponds with the 
unstable fundamental periodic orbit of the Henon-Heiles 
problem near its ground state - in the {^a, ^b) plane the 
projection of an unbalanced orbits is roughly a parabola 
that does not pass through the projection of the AFM 
fixed point. Like the first two orbit families, the behav- 
ior of one spin in the unbalanced orbit is different from 
the other two; therefore there are three unbalanced orbits 
for each of the two AFM fixed points, for a total of six 
unbalanced orbits. At low energies, the odd spin moves 
along a closed curve in {(/>, 6) space while the other two 



spins move along open curves which are dented on one 
corner and are mirror images of one another. One unbal- 
anced orbit lies on the line Zb = line of the Za = 
surface of section. 

The unbalanced orbit disappears at the polar transi- 
tion, at which point the trajectory of the odd spin grazes 
the poles, touching the holes that appear between Fig. 
|c and |d. 



E. Dynamics near the ferromagnetic end 

The periodic orbits are the "skeleton" of the dynamics 
of a system: to understand the "flesh" requires the global 
view obtained through Poincare sections. Choosing a 
good trigger plane for our section was a matter of study- 
ing the projection of orbits in the {^a, ^b) plane: to en- 
sure that all fundamental orbits appear in the Poincare 
section, our criteria were that: (i) all orbits crossed the 
trigger plane, and (ii) no orbits were confined to the trig- 
ger plane, ^a + ^b = satisfied both requirements. 

($, Z) coordinates improve the quality of our Poincare 
sections compared to previous works on the three-spin 
system. Q In previous works, Poincare sections were 
taken with trigger dSiz/dt = and projected on the Sx 
and Sy planes. When this is done, the collective preces- 
sion of the three spins cannot be visually separated from 
more interesting degrees of freedom. Although the con- 
centric loops of KAM tori can be seen in the figures of 
Q when the trajectories on the tori are non-precessing, 
they are superimposed by random dots from precessing 
chaotic trajectories. Worse, at energies close to the an- 
tiferromagnetic ground state {E = —1.5), trajectories on 
KAM tori themselves precess, destroying their image. As 
a result, the sections of jj] had limited utility as a map of 
the dynamics of the three spin system and had to be sup- 
plemented with power spectra of the classical trajectories 
to determine if trajectories were regular or chaotic. 

Fig. |l2| is a surface of section using the $^ + ^b — 
trigger which we produced using a method of visualiz- 
ing Poincare sections for two degree of freedom systems 
in three dimensional space described in Appendix 
Fig. |l^ is a single image of a simulated 3 dimensional 
object which can be interactively rotated and viewed 
from arbitrary positions. The dark opaque object is 
the E = 2.05 energy surface, which is approximately an 
oblate spheroid with the {Za, Zb) plane passing through 
the equator. Over that surface is plotted cloud of dots 
which are the intersections oi E ~ 2.0 trajectories with 
the surface of section. 

All of the fundamental periodic orbits intersect the sur- 
face of section in two places, once passing through the 
surface of section in the positive direction {^a > 0) and 
once in the negative direction {^a > 0.) Just on the 
lower visible edge of the energy surface is a sort of ter- 
minator which divides trajectories that cross the surface 
of section in the positive and negative directions; this 
curve is not quite a geodesic but it does divide the energy 
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surface into two approximate hemispheres. The intersec- 
tions of the stationary-spin and counterbalanced orbits 
with the surface of section form a ring of 12 fixed points 
lying in the {Za,Zb) plane with exact 12-fold symme- 
try while the two three-phase orbits cross the surface of 
section away from the plane. 

The nature of KAM tori in the ferromagnetic limit 
— > 3 is visible in Fig. |2|. A concentric family of KAM 
tori exist around each counterbalanced orbit, and families 
of tori also exist centered around the three-phase orbits. 
Stationary spin orbits lie on the separatrix which divides 
counterbalanced tori from three spin tori. As energy is 
lowered, this separatrix is the first place where tori break 
and chaos is observed. 



F. Dynamics near the antiferromagnetic end 



To study the dynamics of the three-spin system near 
the antiferromagnetic end, i5 ^ — | we chose <i>^ = as 
a tr igger. Although this violates criterion (ii) of Section 
VE, we gain the advantage that this trigger plane ex- 



tends from the lowest to the highest energies and crosses 
both antiferromagnetic fixed points. The practical dis- 
advantage is that a stationary spin orbit exists on the 
= line, and appears as a curve on the Poincare 
section rather than a point, but this does not terribly 
complicate the interpretation of the section. 

Fig. [l^ illustrates the transition to chaos in the anti- 
ferromagnetic regime with a = 0.5. At E = —1.39 (see 
Fig. 14) most tori are unbroken and motion is primarily 
regular. By E — —1.35 chaos is becoming noticeable in 
separatrix regions, and hy E — —1.3 chaos is widespread. 
At E = —1.2 no islands of regular motion are obvious. 
However, we know that regular islands do exist because 
the inner branch of the stationary spin orbit is stable at 
some energies in this regime (See section ^) The tran- 
sition to chaos on the antiferromagnetic side has been 
observed previously Q in the same energy range. 



G. Symmetry and dynamics 



Because the Hamiltonian ( |2.3| ) is threefold symmetric 
arou nd the antiferromagnetic fixed points, (see Section 
II B) the low energy behavior of the three spin system 
falls into the same universality class as the well-known 
Henon-Heiles system with Hamiltonian |l^ 

H = ^{pl+pl)+^{x'+y')+x{^'y~y'm. (5.2) 



This can be seen in Fig. 14, which looks remarkably like 
Poincare sections of the Henon-Heiles system. 

Symmetries around a fixed point determine many 
properties of the dynamics of a system in its vicinity in- 
cluding the nature of the fundamental orbits, the global 
geometry of trajectories in phase space, and degenera- 
cies in orbit frequencies. Threefold rotation symmetry 



around the antiferromagnetic fixed point ensures that the 
periodic orbits and KAM tori near one AFM fixed point 
can be mapped 1-1 to those in Henon-Heiles, but it also 
guarantees that the frequencies of all periodic orbits con- 
verge in the E —^ —1.5 limit: if we perform a Taylor series 
expansion of the Hamiltonian at the fixed point (as in Eq. 
( |3.5[ )) the only quadratic term compatible with threefold 
rotation symmetry is that with circular symmetry in the 
{^a,Sb) and {Za,Zb) planes. 

The three-spin system exhibits a six-fold rotation sym- 
metry near the FM limit which is responsible for a differ- 
ent orbit and torus geometry in the _E — *■ 3 limit which is 
probably generic for Hamiltonian fixed points with six- 
fold symmetry. With six-fold symmetry, the first and 
second derivatives of the time-energy curves are the same 
for all fundamental orbits at i? = 3; as a result, the peri- 
ods of orbits are remarkably degenerate for a large range 
in energy (see Fig. 0.) 



VI. CONCLUSION 

This paper is a detailed analysis of the classical dynam- 
ics of the three spin cluster with Hamiltonian Eq. ( |l.l| ) ; 
many of the features we find are connected with quantum 
phenomena in the accompanying paper ^ . One class of 
phenomena are connected with changes in the topology 
of the energy surface (see Section IV) which occur as a 
function of energy: if we make a plot of quantum energy 
levels as a function of tr, shown in Fig. 5 of ||l| , we observe 
a tunnel splitting as pairs of near-degenerate levels (at 
E < Ec{a)) cross the E^ia) curve (see Eq. (|^.) This 
is caused by tunneling between quantum states localized 
on the two disconnected parts of the energy surface. A 
second phenomenon related to the topology of the en- 
ergy surface is that the classical and quantum densities 
of states are apparently constant as a function of e nergy 
for Ec{a) < E < -I (see Fig. 4 of [Q and Section |lVD 
of this paper.) 

Our analysis of fundamental periodic orbits in Section 
IV also has significance for the quantum problem. The 
Gutzwiller trace formula |l^ predicts that classical pe- 
riodic orbits cause oscillations in the density of states. 
In we observed these oscillations by applying spectral 
analysis to the quantum density of states, (see Fig. 3 of 

§•) 

Some technical aspects of the work presented in this 
paper are interesting. Fi rst, t he change of coordinates 
that presented in Section II A enables us to understand 
the three-spin cluster better than previous studies |^ 
as we take advantage of the clusters conservation 
of total Sz to reduce the dynamics to a tractable two- 
degree of freedom system. Second, our use of three- 
dimensional visualization for visualizing the energy sur- 
face and Poincare sections clarifies interpretation of 
Poincare sections when the topology of the energy surface 
is complicated. Even in situations where the topology 
is simple (such as is discussed in Section |VE|,) three- 
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dimensional visualization hides fewer symmetries of sys- 
tem than the customary two-dimensional projection and 
eliminates the confusion caused when two sheets of the 
energy surface are projected on top of one another. This 
method is described further in Appendix ^ 

The three-spin cluster has similiaries to certain well- 
studied systems. Our system has features in common 
with the well-known Anisotropic Kepler Problem (AKP). 

. Like the AKP, total angular momentum around the 
z-axis ^ Sf is conserved, leaving two nontrivial degrees 
of freedom. Both the AKP and Eq. (IT) have a single 
parameter (cr in the case of our system) and are noninte- 
grable for all values of the parameter save one (cr = in 
our case.) In both systems, all three classical frequencies 
are identical in the integrable case. For the AKP, the 
integrable case is the ancient Kepler problem in which 
all trajectories are closed ellipses. In our problem, in the 
integrable case (u = 0), all three spins precess around 
the total spin vector Si at a rate proportional to the 
length of the total spin vector. Our system is dif- 
ferent from the AKP in a number of ways. First, the 
AKP is highly chaotic throughout the parameter space 
in which is has been studied [|l8j (The first stable peri- 
odic orbit was found after the AKP had been studied for 
14 years. ||2^) Our system, on the other hand, shows 
highly regular behavior in much of the parameter space 
(For instance, when > 2 in the a — 0.5 case) as well as 
irregular behavior in other areas (For instance, tr — 0.5 
and E K, —1.2.) Thus the elegant application of symbolic 
dynamics to the AKP ||l^ is not possible for our system. 

Another connection between the three-spin cluster and 
a well-studied system is the similarity between the dy- 
namics of the three-spin cluster in the antiferromagnetic 
limit and the Henon-Heiles problem. |l^, |l] The connec- 
tion here is most obvious in the surface of section shown 
in Fig. 1^ and occurs because the three-spin cluster has 
a three-fold rotational symmetry around the antiferro- 
magnetic ground states similar to the symmetry of the 
Henon-Heiles problem. 

In this work we have gotten a more intimate under- 
standing of a nonintegrable spin cluster than has been 
previously available supporting the work described in [Q , 
which establishes that periodic orbit theory can be ap- 
plied to spin. This work was funded by NSF Grant DMR- 
9612304, using computer facilities of the Cornell Center 
for Materials Research supported by NSF grant DMR- 
9632275. We would hke to thank Masa Tsuchiya, Jim 
Sethna, and Greg Ezra for interesting discussions. 



APPENDIX A: NUMERICAL INTEGRATION IN 
•l", Z COORDINATES 

An important decision in the numerical study of the 
three spin problem is the choice of variables to used to 
integrate the equation of motion. This choice affects the 
speed, complexity and reliability of integration as well as 
the range of Poincare sections that can be easily taken. 



Our ODE integrator library was written in Java and 
evolved from the software used for the results published 
in For both vector components and {^,Z) coor- 

dinates we used adaptive fifth-order Runge-Kutta inte- 
gration based on the code from |^ although our system 
allows the use of different integrators such as fourth-order 
fixed Runge-Kutta for testing. |2^ 

In the early phase of this work we integrated the spin 
vector components {Sx, Sy, Sz) of the individual spins. 
The vector component representation has several ad- 
vantages: the software requirements are simple and it's 
straightforward to write a general routine for evaluating 
the equations of motion for any spin Hamiltonian which 
is polynomial in Sx, Sy and Sz- Spin vector coordinates 
are also free of obnoxious singularities. However, the 
need to isolate the overall precession of the spins from 
more interesting motions led us to integrate the system 
in ((f>, Z) coordinates so we could easily set Poincare sec- 
tions in the $-space. (the rationale for setting triggers in 
$ space is discussed in Section VE.) 

Although the transformation (0, z) — > ($, Z) is a 
straightforward linear transformation, the need to use in- 
verse trigonometric functions to convert S into (0, z) adds 
overhead and, more seriously, additional complexity to 
deal with branch cuts. (Numerical algorithms that work 
with branch cuts, particularly involving square roots, are 
difficult to design. Failure modes caused by roundoff er- 
ror with a probability of 10~^ per dynamical timescale 
are a major complication for a program that calculates 
thousands of trajectories.) Simple strategies for disam- 
biguating branch cuts that do not introduce an error- 
prone memory between steps lose valuable topological 
information. If, for instance, $o is computed from some 
manipulation of the S,; components, and, say, is always 
in the range < $o < 47r\/3 it isn't as easy to determine 
the net precession (secular trend of $0) of a periodic orbit 
as it would be if $0 were integrated directly. 

We performed all of the integrations in this work in 
(<&, Z) coo rdin ates, with equations of motion derived 
from Eq. (2_^). The main difficulty we had is that the 
integration can fail on a trajectory on which a spin passes 
through a pole; this is not a failing of the (<&, Z) coor- 
dinates as much as of the (0, z) coordinates. As spin 
i passes close to the pole, the singularity in the map- 
ping from Zi) Si forces (pi 00. If the trajectory 
misses the pole by more than le — 6 in the z axis when 
using double precision math, the primary consequence is 
that the adaptive step size integrator reduces the time 
step and integration is slowed. If the trajectory passes 
much closer to the pole, however, no step size may have 
a sufficiently small error estimate and the adaptive step 
size algorithm will fail. Although it would be possible 
to avoid this problem by either switching to spin vec- 
tor coordinates when the trajectory passes close to the 
pole or by adding more intelligence (and possibly bugs) 
to the adaptive step size algorithm, in practice it affects 
a small enough volume of phase space that it only man- 
ifests when investigating trajectories specifically chosen 
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to pass near a pole. 

Sometimes it is necessary to work with nor- 
malized to the unit hexagon, for instance, to set a 
Poincare section trigger on $^ = 0. At energies above the 
upper polar transition and below the lower polar tran- 
sition (see Section IV) this is not necessary because the 
trajectory does not wander long distances in the {^a,^b) 
plane. For some time this prevented us from taking 
Poincare sections in the region between the two tran- 
sitions, since the trajectory would eventually wander far 
from the trigger plane. To solve this, we found an algo- 
rithm for mapping (^a,^b) back to the unit hexagon: 
first, (i) use the modulus function to map points into the 
primitive cell of the hexagonal lattice, a rhombus. Then, 
(ii) apply a unit vector translation to those points that 
fall on corners of the rhombus outside the unit hex to 
bring them into the unit hex. 



APPENDIX B: SOLID POINCARE SECTIONS 

In the process of studying the three spin problem, 
we found conventional methods of drawing Poincare sec- 
tions inadequate and improved upon them by develop- 
ing a method for rendering Poincare sections in three- 
dimensional space. This greatly simplified the interpre- 
tation of Poincare sections for our system. Although not 
all systems pose as serious technical problems as ours, 
we believe that this method clarifies the geometry of 
Poincare sections and can simplify the presentation of 
Poincare sections to audiences which are not specialized 
in dynamics. The techniques described in this appendix 
were used to generate Fig. ^, Fig. |l2| and Fig. |l^. 

For Hamiltonian systems with two degrees of freedom, 
the intersection of the energy surface with the surface 
of section is a two-dimensional surface embedded in a 
three-dimensional space (the surface of section). Often, 
this intersection has the topology of a sphere - this is 
true of the Henon-Heiles system as well as for the three 
spin systems above the upper polar threshold and below 
the coalescence energy Ec{<t). The trajectory crosses the 
surface of section in two directions, which we will call 
the positive and negative directions. Over part of the 
sphere, the trajectory crosses in the positive direction 
and over the rest of the sphere, the trajectory crosses in 
the negative direction. In between there is a seam over 
which the trajectory is tangent to the surface of section. 

Difficulties arise when plotting the intersection of tra- 
jectories with the surface of section even when the 
topology is simple. To take a specific example, con- 
sider the case of the three spin problem with a = 0.5 
and E = —2 with the trigger on $^ -I- <i>s = 0. A 
reasonable set of {x,y,z) variables for this problem is 
(x, y, z) — {Za, Zb, ^a)- Conventional choices of projec- 
tion are (x, y) which would overlap the two distinct three 
phase orbits and place the stationary spin orbits on the 
edge of the plot, or (cc, z) which obscures the symmetry 
of the stationary spin and counterbalanced orbits seen in 



the (a;, y) plane. 

While writing the paper we began development 
of a set of Java class libraries for integrating differential 
equations and taking Poincare sections. Interchangeable 
trigger modules are connected to the differential equa- 
tion integrator, allowing the user to set a triggering cri- 
teria for surface of section. When the surface is crossed, 
we solve for the crossing time by varying the integration 
step St and solving for the value of 6t that intersects the 
section by the Newton-Raphson method. (We also tried 
the method of method of Henon ||2^ in which a change 
of variable is made to make the trigger coordinate be- 
come a time coordinate - we found that performance of 
our method and Henon's method is similar, but ours was 
more robust) We primarily integrate with an adaptive 
fifth-order Runge-Kutta implementation, but the design 
of our program allows us to replace our integrator with 
another fixed or adaptive step algorithm. 

One component of a 3-d Poincare section is a rendering 
of the energy surface. For an arbitrary system, the en- 
ergy surface can be rendered by treating it as an implicit 
function; 



E{x,y,z) = Eq 



(Bl) 



this sort of surface can be rendered using ray-tracing 
techniques, or can be rendered in a primitive manner by 
dividing the 3-space into voxels and coloring in cubes 
which are above or below a threshold value. A better 
method is to generate a polygonal mesh with Bloomen- 
thal's algorithm ||2^ . The relative merits of various meth- 
ods for polygonizing a mesh are discussed in j2^. For a 
uniform mesh (not adaptively sampled), Bloomenthal's 
algorithm is equivalent to a table-driven method known 
as Marching Cubes p7| . We can compute a mesh for our 
system in under 20 seconds and 2 megabytes of storage, 
and this time could be reduced by the use of Marching 
Cubes. The resulting mesh can be rendered with a stan- 
dard 3-d system, either an interactive system such as as 
OpenGL or VRML or an off-line system such as POVRay 
or Renderman. 

The generation of the actual Poincare sections is quite 
straightforward. We merely plot a set of (x,y,z) points 
when the trajectory intersects the surface of section. The 
last problem is determining initial conditions for injec- 
tion. In some cases we wish to interactively choose in- 
jection points, in which case it's necessary to translate 
a hand-chosen injection point in screen cordinates into 
into a point one on the energy surface; this way, in an 
interactive 3-d environment, a user can click on the en- 
ergy surface to set an injection point. For most cases, 
we prefer to have the computer automatically generate 
Poincare sections, choosing injection points that reveal 
all major phase space structures. To accomplish this, our 
program divides the surface of section into voxels (tiny 
cubes) and inspects the value of the energy at the corners 
of the cubes to find which cubes intersect the energy sur- 
face. Next, the program scans the surface, ensuring at 
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least one trajectory either originates in or passes through 
each voxel that contains part of the energy surface. To 
inject in a voxel, our program computes the gradient of 
the Hamiltonian at the center of the voxel and precisely 
locates the energy surface by a Newton search along the 
line of fastest change in energy. It then marks that voxel 
as filled. The trajectory is then evolved forward in time 
for a fixed number of intersections: at each intersection, 
the containing voxel is set as filled. The program con- 
tinues to find voxels which are not filled and injects into 
them until all voxels are filled. 

Although the energy surfaces are interesting in them- 
selves for systems with topology changes (such as the 
three spin system), they are also essential for interpret- 
ing solid Poincare sections. To make a solid Poincare sec- 
tion we render a cloud of dots in three-dimensional space 
at the points where trajectories intersect the surface of 
section. The resulting dot cloud is transparent and is 
difficult to make sense of without an opaque object un- 
derneath it. Several factors conspire to confuse the eye. 
For instance, a part of a dot cloud that is further from 
the viewer is made visually denser by perspective which 
makes it appear heavier and more solid to the viewer - 
causing the viewer to conclude that the distant part of 
the dot cloud is closer. An opaque rendering of the en- 
ergy surface eliminates the confusion caused when images 
of the energy surface overlap. 

Because integrating differential equations is slow, it 
takes about two hours to compute our Poincare sections. 
Because the calculation can be broken into a a number 
of mainly unrelated calculations, parallelizing the calcu- 
lation across multiple processors on an SMP machine or a 
cluster of computers would be straightforward and worth- 
while. 

A complicated aspect of the three-spin problem is the 



hexagonal tiling of the phase space. This introduced two 
difficulties which required special solutions. We imple- 
mented these solutions using object-oriented techniques 
to replace general-purpose classes with specialized sub- 
classes. First, the 3-space {Za, Zb,^a) is bounded by 
the points at which the spins touch the north pole. The 
volume insid e this boundary is a hexagonal prism (See 
Section II C and Fig. ^) Using a cubical grid to rep- 
resent a non-cubical domain forced us to keep track of 
which points were inside the hexagonal prism and which 
were outside. Not only was this difficult and slow, but 
the irregular pixilation of the boundary confused both 
the surface-constructing algorithms and the eye of the 
viewer, resulting in surfaces with a jagged, serrated ap- 
pearance. We solved the problem by replacing the Java 
class responsible for mapping integer voxel numbers to 
(x, y, z) coordinates with one that smoothly maps a cubic 
grid on an integer lattice to a hexagonal prism. All sam- 
ple points fell within the allowed range. We headed off 
numerical problems caused by square roots in the energy 
function by multiplying x and y coordinates by 1 — 10~^^ 
before computing the energy. 

At energies between the upper and lower polar thresh- 
olds, another set of problems arises. Here, the projection 
of a trajectory in {^a, ^b) can wander outside the unit 
hexagon. If the trajectory wanders in the {^a, ^b) plane 
it may fail to intersect the surface of section. This prob- 
lem can be solved by wrapping the Java class which pro- 
vides for inspection and manipulation of the system state 
in ($, Z) coordinates (the actual integrator was written 
before these coordinates were settled) with a Proxy 
class which maps {^,Z) coordinates back to the unit 
hexagon - an algorithm for performing this mapping is 
described in Appendix |^. 
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rows are the three lattice vectors. The primitive unit cell is 
solid and copies are dashed. The following fixed points are 
depicted: (i) the one ferromagnetic (FM) fixed point, (ii) the 
two antiferromagnetic (AFM) fixed points, and (iii) the three 
antiparallel (A) fixed points. 
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FIG. 3: Energies at which transitions in the energy surface 
topology occur. 
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FIG. 4: Changes in the topology of the energy surface 
{H{Za,Zb,^b) = E with $a = 0), for a > 2/3, here 
a = 0.75. (a)£ = -0.92, (b) £ = - = 0.98, (c) E = -1.02, 
and (d) E = —1.1. The hexagonal prism represents the hexag- 
onal boundaries of the {Za,Zb) plane: the horizontal direc- 
tion is Za , the vertical direction is Zb and the long axis of the 
prism is $s. The endpoints of the prism are at $_b = ^^vr 
and $B = i^TT. 





FIG. 5: The trajectory of a moving spin on the stationary 
spin manifold, for a = 0.5. The well around the negative z- 
axis is the antiparallel fixed point, while the wells above and 
below it are the two antiferromagnetic fixed points. (The FM 
fixed point is on the positive x-axis and is on the hidden side 
of the sphere.) The separatrix between orbits centered on the 
AFM fixed points those centered aro und the antiparallel fixed 
points occurs at energy Ec (see Eq. (4.2).) 
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FIG. 6: The classical density of states as a ftmction of energy 
for a = 0.6 determined by numerical integration, (a) the 
coalescence transition, (b) the lower polar transition, (c) the 
antiparallel transition. The total area under the curve is Stt, 
the volume of the reduced phase space. 
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FIG. 7: Energy-period curve of the three spin system with 
a = 0.5. (a) Stationary Spin, (b) Counterbalanced, (c) Three- 
phase orbit, and (d) Unbalanced orbits. 
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FIG. 8: The stability exponent p (p = Tr M where M is the 
stabiUty matrix) of the stationary spin orbit as a function of 

energy for a = 0.5. An orbit is stable when \p\ < 2 and unsta- 
ble when \p\ > 2. The stationary spin orbit is stable near the 
AFM ground state and alternates between being stable and 
unstable until Ec- The outer branch of the orbit is unstable 
for all E > Ec while the inner branch appears to always be 
stable until its disappearance at E = 1. 
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FIG. 9: The stability exponent p = TrM where M is the 
stability matrix for the counterbalanced orbit in the case a = 
0.5 where the orbit is always stable. 



16 



30.0 

in 



2 20.0 - 



Q. 

10.0 




-1.50 -1.25 -1.00 -0.75 

energy 

FIG. 10: The precession rate squared as a function of energy 
for the three phase orbit with cr = i. The precession rate is 
the amount that $o increases after one period of the orbit in 
the reduced space. In the Hmit of i? — > — | the precession rate 
becomes 2n. The straight line fit (fit to all of the points with 

E > 0.90) indicates {E - Et,)^ scaUng with Eb = -0.7547. 
Note that {2-!vf « 39.48. 




FIG. 11: Stability exponent p for the three phase orbits for 
a = 0.5. The three phase orbit is stable for E > Eb, Eb ~ 
0.75. Starting from high energies, p reaches 2 at Eb where 
a bifurcation gives birth to two stable precessing orbits and 
an unstable nonprecessing orbit. Another (unstudied) bifur- 
cation happens around E = 0.91. Starting at the low-energy 
AFM end, the orbit starts out stable and becomes unstable 
at E = —1.Z2 in a period doubling bifurcation. 



FIG. 12: Poincare section ai a — 0.5, E — 2.0 with a $a + 
~ trigger: point (a) is a stationary spin orbit, point (b) 
is a counterbalanced orbit, and point (c) is a three-phase orbit. 
The energy surface is approximately an oblate spheroid with 
stationary spin and counterbalanced orbits on the equator on 
which <1>B =0; the back side of the surface is hidden. Just 
behind the equator is the other three-spin orbit for which the 
spins rotate opposite the direction of the visible three-spin 
orbit. Throughout most of the visible hemisphere, orbits cross 
the section traveling in the positive direction (t&A > 0); on 
the invisible hemisphere, the picture is mirror-reflected across 
the "terminator" for orbits crossing in the negative direction. 
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FIG. 13: The transition to chaos at the antiferromagnetic end 
[E —1.5) with a — 0.5. Shown is one of the two lobes from 
Fig. ^. See Fig. |l^ to determine which orbits are which. 




FIG. 14: A projection of the E = —1.39 Poincare section 
into the ($b,Zb) plane. Orbits marked (a) are three phase 
orbits, (b) is a stationary spin orbit, and (c) is an unbalanced 
orbit. 



